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Q\ | Abstract 

The main concern of the present paper is to study large-time behavior of 
solutions to an ideal polytropic model of compressible viscous gases in one- 
, dimensional half space. We consider an outflow problem, where the gas blows 

7— I ' out through the boundary, and obtain a convergence rate of solutions toward 

a corresponding stationary solution. Here the existence of the stationary 
solution is proved under a smallness condition on the boundary data with the 
. ^ | aid of center manifold theory. We also show the time asymptotic stability of 

^ ■ the stationary solution under smallness assumptions on the boundary data and 

the initial perturbation in the Sobolev space, by employing an energy method. 
Moreover, the convergence rate of the solution toward the stationary solution 
is obtained, provided that the initial perturbation belongs to the weighted 
Sobolev space. Precisely, the convergence rate we obtain coincides with the 
spatial decay rate of the initial perturbation. The proof is mainly based on a 
priori estimates of the perturbation from the stationary solution, which are 
derived by a time and space weighted energy method. 
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1 Introduction and main result 
1.1 Formulation of the problem 

We study large-time behavior of a solution to an initial boundary value problem for 
the compressible Navier-Stokes equations over one-dimensional half space M + := 
(0, oo). An ideal polytropic model of compressible viscous fluid is formulated in the 
Eulerian coordinates as 

pt + (pu) x = 0, (1.1a) 
(pu) t + {pu 2 + p(p,9)) x = pu xx , (1.1b) 

{p(c v 6 + y)} t + {p u ( c vd + y) +p(p> 9 ) u } x = (puu x + n9 x ) x , (1.1c) 

where unknown functions are p = p(t,x), u = u(t,x) and 9 = 9(t,x) standing for a 
mass density, a fluid velocity and an absolute temperature, respectively. Due to the 
Boyle-Charles law, a pressure p is explicitly given by a function of the density and 
the absolute temperature: 

p = p(p,9):=Rp9, 

where R > is a gas constant. Positive constants c v , p and k mean a specific heat 
at constant volume, a viscosity coefficient and a thermal conductivity, respectively. 
Due to Mayler's relation for the ideal gas, the specific heat c v is expressed by the 
gas constant R and an adiabatic constant 7 > 1 as 

R 

7-1 

We also introduce physical constants 

Cp . 7^v T^i -^r ' T^"> 

7 — 1 K K7 - 1 

which stand for a specific heat at constant pressure and the Prandtl number, re- 
spectively. The Prandtl number plays an important role in analysis of a property of 
a stationary solution. 

We put an initial condition 

(p,u,6)(0,x) = ( Po ,u ,6 )(x) (1.2) 

and boundary conditions 

u(t,0)=u h <0, 9(t,0) = 9 h >0, (1.3) 

where Mb and 6\> are constants. It is assumed that the initial data converges to a 
constant as x tends to infinity: 

lim (p ,u ,9o)(x) = (p+,u + ,9+). 

x— >oo 

Moreover, we assume that the initial density and absolute temperature are uniformly 
positive, that is, 

inf p (x) > 0, inf 9 (x) > 0, p + > 0, 9+ > 0. 
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The boundary condition for u in (11.31) means that the fluid blows out from the 
boundary. Hence this problem is called an outflow problem (see [H]). Due to the 
outflow boundary condition, the characteristic of the hyperbolic equation (j 1.1 aft for 
the density p is negative around the boundary so that two boundary conditions are 
necessary and sufficient for the wellposedness of this problem. 

In the paper [9], Kawashima, Nishibata and Zhu considered the outflow problem 
for an isentropic model and obtained a necessary and sufficient condition for the 
existence of the stationary solution. Moreover, they proved the asymptotic stability 
of the stationary solution under the smallness assumption on the initial perturbation 
and the strength of the boundary data. A convergence rate toward the stationary 
solution for this model was obtained by Nakamura, Nishibata and Yuge in [TJ] 
under the assumption that the initial perturbation belongs to the suitably weighted 
Sobolev space. The main concern of the present paper is to extend these results 
to the model of heat-conductive viscous gas. Precisely, we show the existence and 
the asymptotic stability of the stationary solution as well as the convergence rate 
for the ideal polytropic model ( 11. ip . Compared to the isentropic model, the heat- 
conductive model is more difficult to handle. For example, since the model (II. ip 
has two parabolic equations, the equations for the stationary wave are deduced to 
a 2 x 2 system of autonomous ordinary differential equations. However, it becomes 
a scalar equation in the case of the isentropic flow. Therefore, to obtain a condition 
which guarantees the existence of the stationary solution for the heat-conductive 
model, we have to examine dynamics around an equilibrium of the system by using 
center manifold theory. 



1.2 Dimensionless form 

For the stability analysis on the equations (II. ip . it is convenient to reformulate 
the problem into that in the dimensionless form. For this purpose, we define new 
variables x and t by 

x - t 

x '- = v t:= r 

where L and T are positive constants. We also employ new unknown functions 
(p, u, 6) defined by 

1 .1 - - 1 

p(t,x) := — p(t,x), u(t, x) := ru(t,x), 9{t,x) := —9(t,x). (1.4) 

P+ ' \u+\ 0+ 

Here we note that the constant u + must satisfy 

m+<0 (1.5) 

for the existence of the stationary solution. Indeed, the stationary solution (p, it, 9)(x) 
satisfies 

p(x)u(x) = p+u + , (1.6) 

which is obtained by integrating (pu) x = over (x, oo). Substituting x = in (11.61) . 
we get u + = p(0)«b/ p+, which immediately yields (I1.5P by using the positivity of the 
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density and the boundary condition u b < 0. Next we define dimensionless physical 
constants by 

P-= — ; «:= — ; 14, c v := — (1.7) 

P+\ u +\ P+\ u +\ 7(7-1) 

and a dimensionless pressure by 

p=p(p, 9) := -p9. 

7 

We also introduce Mach number M + at the spatial asymptotic state: 

NU := H 



where c + := ^R^9 + is sound speed. Using the dimensionless constants (jl.7p . we 
represent the Prandtl number P r as 



-Pr " 



kMKj- 1)' 



Substituting (11. 4p in (11.11) and letting L = and T = 1, we have the equations 
for (p, u, 9) in the dimensionless form as 

pt + (pu) x = 0, (1.8a) 
Wi + [pu 2 + JpP(P> d )) x = P u ^ ( L8b ) 

(1.8c) 

In the equations (jl.8p . without any confusion, we abbreviate the symbol " to 
express dimensionless quantities. The initial and the boundary conditions for the 
dimensionless function (p, u, 9) are prescribed as 

Po _uo_ Oo 

,p+ \u+y # H 

\im(p ,u ,9 )(x) = (1,-1,1), (1.9b) 



(p, u, 0)(O, x) = (p , u , 9 )(x) := f — , 7— —7 j ^) (ar), (1.9a) 



(„ > fi)(t,0) = (« b A):=(A,r-)- (1.10) 

We also abbreviate the hat and write the dimensionless initial data and boundary 
data as (po,uo,9 Q ) and (w b , #b) respectively in (II. 9p and (11.101) . 



1.3 Main results 



The main concern of the present paper is to consider the large-time behavior of 
solutions to the problem (11.81) . (II. 9p and (jl.lOp . Precisely we show that the solution 
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converges to a stationary solution (p,u,9)(x), which is a solution to (11.81) indepen- 
dent of time variable t. Thus the stationary solution (p, u, 9) satisfies the system 

(pu)x = 0, (1.11a) 
PU 2+ M2P) X = l™**> (1.11b) 

{fa ( a^ 6 ' + y ) + ~w^}x = + K ^ a " (l.iic) 

where p := p(p, 0). The stationary solution is supposed to satisfy the same boundary 
condition (11.101) and the same spatial asymptotic condition ( 11.9bl) : 

(u,9)(0) = (iibA), lim(p,M)(x) = (1,-1,1). (1.12) 

We summarize the existence and the decay property of the stationary solution 
(p, u, 9) satisfying (II. lip and (I1.12p in the following proposition. To this end, we 
define a boundary strength 5 as 

6:= |K + 1A-1)|. 

Proposition 1.1. Suppose that the boundary data (wbj^b) satisfies 

(u h ,9 h )eM + :={(u,9)em 2 ; |(«+l,0-l)|<e o } (1.13) 

for a certain positive constant e . Notice that the condition (11.131) is equivalent to 
5 < e . 

(i) For the supersonic case M + > 1, there exists a unique smooth solution (p,u,9) 
to the problem (II. lip and f 1 1 . 1 2 j) satisfying 

\d*(p(x) - l,u(x) + l,9(x) - 1)| < C5e~ cx for k = 0, 1, 2, . . . , (1.14) 

where C and c are positive constants. 

(ii) For the transonic case M + = 1, there exists a certain region Ai° C M. + such 
that if the boundary data («b, 9b) satisfies the condition 

(u h ,9 h )eM°, (1.15) 

then there exists a unique smooth solution (p, u, 9) satisfying 

\d k x {p{x) - 1, u(x) + 1, 9(x) - 1)| < C ^\ k+1 + C5e~™ for fc=0,l,2,.... 

(1.16) 

(iii) For the subsonic case M + < 1, there exists a certain curve M.~ C M. + such 
that if the boundary data (u^, 9^) satisfies the condition 

(u h ,9 h )eM-, (1.17) 

then there exists a unique smooth solution (p,u,9) satisfying (11.141) . 
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Figure 1: For the transonic case M + = 1, the region A4° consists of one side of A4 + 
divided by the local stable manifold 9 = h s (u). For the subsonic case M + < 1, the 
curve M.~ coincides with the local stable manifold. 

The rough sketches of the regions Ai + , A4° and M.~ are drawn in Figured! The 
precise definitions of A4° and Ai~ are given in (I2.19p . The boundary of A4°, which 
is the stable manifold for the stationary problem, is a curve in the state space. The 
geometric property of this curve is completely characterized by the Prandtl number. 
This observation is discussed in Section 12.31 

The asymptotic stability of the stationary solution (p, u, 6) is stated in the next 
theorem. 

Theorem 1.2. Suppose that the stationary solution (p,u,6) exists. Namely it is 
assumed that one of the following three conditions holds: (i) M + > 1 and ( 11.131) , 
(ii) M + = 1 and ( 11.151) . (iii) M + < 1 and (11.171) . In addition, the initial data 
(po, ^o, #o) is supposed to satisfy 

( Po ,u ,9 )-(p,u,9)eH 1 (R + ) 
for a certain constant a G (0, 1). Then there exists a positive constant E\ such that 

||(po,MoA) - (p,u,9)\\ H i +5 < ei, 
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then the initial boundary value problem ( ll.8p . ( I1.9P and (11.101) has a unique solution 
globally in time satisfying 

(p, u, 0) - 05, u, 9) E C([0, oo); ^ 

for an arbitrary T > 0. Moreover, the solution (p,u,8) converges to the stationary 
solution (p, u, 9) uniformly as time tends to infinity: 

lim ||0, - {p,u,6)\\ L °° = 0. (1.19) 

t— >oo 

We also show a convergence rate for the stability (11.191) by assuming additionally 
that the initial perturbation belongs to the weighted Sobolev space. 

Theorem 1.3. Suppose that the same conditions as in Theorem \1.2\ hold. 

(i) For the supersonic case M + > 1, if the initial perturbation satisfies 



(po,u o ,0 o ) - (p,u, 9) E L 



2 , 



for a certain positive constant a, then the solution (p,u,8) to ( 11. 81) . ( 11. 91) and 
( 11.101) satisfies the decay estimate 

i,u,e){t) - (p,u,e)\\L°° < c(i + i)- Q / 2 . (i.2o) 



(ii) For the transonic case M + = 1, let a E [1,2(1 + \/2)). There exists a positive 
constant e 2 such that if 

6- 1 / 2 \\(p ,u ,9 )-(p,u,9)\\ Hi <e 2 , 

then the solution (p, u, 9) satisfies the decay estimate 

\\{p,u,0){t)-{p,uM\L~ <C{l + t)- a '\ (1.21) 

Remark 1.4. (i) For the supersonic case M + > 1, we can prove an exponential 
convergence rate 

\\(p,u,9)(t)-(p,u,9)\\ Lao <Ce- at 
provided that the initial data satisfies the conditions as in Theorem 11.21 and 

( Po ,«oA) - (~P,u,9) E Lj?, exp (R+) := {u E Ll c (R + );eW*u E L 2 (R + )}, 

where a is a positive constant depending on (. Since the proof is almost same as 
that for the isentropic model studied in the paper [13], we omit the details. 

(ii) To obtain the convergence rates (11.201) and (11.211) . we derive weighted energy 
estimates. In the derivation, we essentially use a property that all of characteristics 
of a hyperbolic system, which is obtained by letting p = and k = in ( 11.81) . are 
non-positive at spatial asymptotic state. However, for the subsonic case M + < 1, 
one characteristic is positive. Due to this, it is difficult to obtain a convergence rate 
for the subsonic case by using the weighted energy method. 

(iii) Compared with the results in [SI H21 US] considering the convergence rate for 
a scalar viscous conservation law, the convergence rates in (ll.20p and (ll.2ip seem 
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optimal. For the transonic case, owing to the degenerate property of the stationary 
solution, the weight exponent a needs to be less than a certain constant, i.e., a < 
2(1 + \/2). This kind of restriction on the weight exponent is also necessary to 
obtain a convergence rate 0(t~ Q / 4 ) toward the degenerate nonlinear waves for a 
scalar viscous conservation law and an isentropic model studied in the papers p2[ 
HH [HI US]- We note that, in the papers [T8J [19], the same restriction a < 2(1 + 
\/2) is also required for an isentropic model and a scalar viscous conservation law 
Ut + f(u) x = u xx , where a degeneracy exponent is equal to 1, that is, f{u) = 
C{u — u + ) 2 + 0{\u — m+| 3 ). Recently, Kawashima and Kurata in [7] studied the 
stability of the degenerate stationary solution for a viscous conservation law and 
obtained the same convergence rate 0(r Q/4 ) by using the weighted energy method 
combined with the Hardy type inequality under a more moderate restriction a < 5, 
which is best possible in the sense that the linearized operator around the degenerate 
stationary solution is not dissipative in L 2 a for a > 5. 

Related results. From the pioneering work (5] by Il'in and Olemik, there have 
been many studies on the stability of several nonlinear waves for a scalar viscous 
conservation law. For instance, Kawashima, Matsumura and Nishihara in P, d2l 
[T5] obtained a convergence rate toward a traveling wave for the Cauchy problem. 
For a one-dimensional half space problem, Liu, Matsumura and Nishihara in [TO] 
considered the stability of the stationary solution. 

For the half space problem of the isentropic model, Kawashima, Nishibata and 
Zhu [9] proved the existence and the asymptotic stability of the stationary solu- 
tion for the outflow problem. The convergence rate for this stability result was 
obtained by Nakamura, Nishibata and Yuge in |JJ:J by assuming that the initial per- 
turbation decays in a spatial direction. The generalization of this one-dimensional 
outflow problem to the multi-dimensional half space problem were studied by Kagei, 
Kawashima, Nakamura and Nishibata in [Bl [T3] . Precisely, Kagei and Kawashima 
in [6] proved the asymptotic stability of a planar stationary solution in a suitable 
Sobolev space. The convergence rate was obtained by Nakamura and Nishibata in 
[T3] . There are also several works on the stationary problem for the Boltzmann 
equation (or BGK model) in half space. See [TJ [2] for numerical computations and 
[T7] for asymptotic analysis. 

Outline of the paper. The remainder of the present paper is organized as fol- 
lows. In Section [2[ we discuss the existence of the stationary solution and present 
the proof of Proposition 11.11 In Section 12. 2\ we show a precise decay property of 
the degenerate stationary solution, which is utilized in the stability analysis of the 
degenerate stationary solution. In Section [3j Theorem 11.21 is proved by deriving 
uniform a priori estimates of the perturbation from the stationary solution in H 1 
Sobolev space by an energy method. Finally, in Section HJ we prove Theorem 11.31 
The crucial argument is to derive time and space weighted energy estimates. For 
the supersonic case, in Section 14.11 we obtain the weighted estimates in L 2 space 
and combine it with the uniform estimates in H l obtained in Section [3j Then we 
obtain the convergence rate (II .201) with the aid of induction. However, owing to the 
degenerate property of the transonic flow, we have to derive the weighted estimate 
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not only in L 2 but also in H 1 in order to obtain the convergence rate (jl.2ip . This 
is discussed in Section 14.21 

Notations. The Gaussian bracket [x] denotes the greatest integer which does not 
exceed x. For p G [l,oo], L P (1R + ) denotes the standard Lebesgue space over M + 
equipped with the norm || • \\lp- We use the notation || • || := || • \\ L 2. For a non- 
negative integer s, H S (M. + ) denotes the s-th order Sobolev space over M + in the L 2 
sense with the norm 



fc=0 

For constants p G [1, oo) and a6l, I^,(R + ) denotes the algebraically weighted 
L p space defined by L£,(R + ) := {u G Lf oc (IR + ) ; \\u\\ l p < oo} equipped with the 



We also use the notation | • \ a := || • \\ L 2. The space H^{R + ) denotes the alge- 
braically weighted H s space corresponding to L 2 a (M. + ) defined by H^(R + ) := {u G 
L 2 t (R + ) ; d k x u G L^(M + ) for k — 0, ... , s}, equipped with the norm 



For a G (0, 1), B a {R + ) denotes the space of the Holder continuous functions over 
M + with the Holder exponent a with respect to x. For a non-negative integer k, 
B k+a (R + ) denotes the space of functions satisfying d l x u G B a (M. + ) for an arbitrary 

1 = 0,...,k equipped with the norm || ■ Hgfc+a. For a, (3 G (0,1) and T > 0, 
£> a,/3 ([0, T] xl + ) denotes the space of the Holder continuous functions over [0, T]xR + 
with the Holder exponents a and (3 with respect to t and x, respectively. For non- 
negative integers k and £, B^ +a ' e+f3 := B k+a,i+,3 ([0, T] x R + ) denotes the space of 
functions satisfying d\u, d x u G B a,,3 ([0,T} x R + ) for arbitrary i = 0, . . . , k and 
j = 0, . . . , £ equipped with the norm || • \\ B k+ a ,e+/3. 

2 Existence of stationary solution 

This section is devoted to showing Proposition 11.11 Precisely we prove the existence 
of a solution to the stationary problem (II. lip and fll,12p . To this end, we reformulate 
the problem (11.1 ID and (I1.12p into a 2 x 2 autonomous system of ordinary differential 
equations of first order. 

2.1 Reformulation of stationary problem 

Integrating (11.11 al) over (x, oo), we have 




norm 





p(x)u(x) 



1. 



(2.1) 
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Integrating fll.llbft and (j 1.1 left over (x, oo) and substituting (12 .ip in the resultant, 
we obtain the system of equations for (u, 6){x) := (u, 9){x) — (—1,1) as 



where J is the Jacobian matrix at an equilibrium point (0, 0) defined by 

J :-- 



-( 1 - -n —L- 

1 _ Cy 



and / and g are nonlinear terms defined by 

t(- a\ u{u + 6) - u 2 

f{u Q) ■= — — , g(u, 9):=—. 

u,M+ r y{u — 1) 2k 

Boundary conditions for (u, 9) are derived from (11.121) as 

(u,9)(0) = (u h + l,9 h -l), \im (u,9)(x) = (0,0). (2.3) 

To prove the existence of the stationary solution (p,u,9), it suffices to show the 
existence of the solution (u, 9) to the boundary value problem (12. 2p and (12. 3p . To this 
end, we diagonalize the system (12.21) . Let Ai and A2 be eigenvalues of the Jacobian 
matrix J. Since we see later that J has real eigenvalues, we assume Ai > A2. 
Let ri and r 2 be eigenvectors of J corresponding to Ai and A2, respectively, and 
let P := (r 1; r 2 ) be a matrix. Furthermore, using the matrix P, we employ new 
unknown functions U(x) and 0(x) defined by 

u ( x )\ ._ p-i fu(x)\ _ , 2A s 



O(x) ) ' \9(x 

We also define a corresponding boundary data and nonlinear terms by 

U b \ _ p -i(u h + l\ ff(U,G)\ _ p -i(f(u,&) 
O h J- ^ \0b-lj' \g{U,e)J- * \g(u,9) 

Using these notations, we rewrite the problem (12.21) and (12. 3p in a diagonal form as 

ig)-ft3(S)+cra. 

(U, G)(0) = (U h , 6> b ), lim (U, 0)(x) = (0, 0). (2.6) 

x— >oo 

Since the existence of the solution to the problem (ll.lip and (11 . 12j) follows from 
that to the problem (12. 5ft and (I2.6p . here we show the latter. Firstly, we consider 
the case M + > 1. Since a discriminant of an eigen-equation of the matrix J satisfies 

(Tr J) 2 - 4 det J = (b - cf + a 2 + 2ab + 2ca > 0, 

where a, b and c are constants defined by 
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the eigenvalues Ai and A2 are real numbers. Moreover we see 

\ l + A 2 = Tr J = -(a + b + c) < 0, AiA 2 = det J = bc> 0, 

which show that Ai < and A 2 < 0. Thus, the equilibrium point (0, 0) of (I2.5P is 
asymptotically stable. Consequently, if |(£/t,j6?b)| is sufficiently small, the problem 
( 12. 5p and (12. 6p has a unique smooth solution (U, 0) satisfying 

\d k x {U{x),0{x))\ < C5e~ cx for k = 0,1, (2.7) 

Next we study the case M + = 1. Since the matrix J satisfies 

c d 

TrJ = — — < 0, detJ = 0, d := u, + 74(7 - l) 2 , 

/IK 

the eigenvalues of J are Ai = and A 2 = —c v d/(u,K) of which eigenvectors are 
explicitly given by 

respectively. Notice that the matrix P = (r 1 ,r 2 ) satisfies det P = —d < 0. Thus 
there exist a local center manifold = h c (U) and a local stable manifold U = h s (0) 
corresponding to the eigenvalues Ai = and A 2 = —c v d/(fiK,), respectively. In order 
to show the existence of the solution, we have to examine dynamics on the center 
manifold. To this end, we employ a solution z = z(x) to (12.51) restricted on the 
center manifold satisfying the equation 

z x = f{z,h c {z)). (2.8) 

By virtue of the center manifold theory in [3], there exists a solution z to (12.81) such 
that the solution (U, 0) to (12.51) and (12.61) is given by 

U(x) = z(x) + 0{5e- cx ), (2.9) 
0{x) = h c {z{x)) + 0(5e~ cx ). (2.10) 

Therefore, to obtain the solution (U,0) to (12.51) and (12. 6p . it suffices to show the 
existence of the solution to (12. 8p satisfying z{x) —> as x — > 00. We see that the 
nonlinear terms / and g satisfy 

f(U, 0) = -^U 2 + O {\U\ 3 + \U0\ + \0\ 2 ) , (2.11) 

9(U, 0) = ^^(P r - 2)U 2 + O (\U\ 3 + \U0\ + \0\ 2 ) . (2.12) 
2\id v ' 

Substituting (I2TTD in (EE}, we deduce (EHD to 

' " 1 z 2 + 0(\z\ 3 ), (2.13) 



2d 

which yields that z is monotonically decreasing for sufficiently small z. Thus, to 
satisfy z(x) — > as x — >• 00, the boundary data z(0) should be positive. Namely, for 
the existence of the solution {U,0), the boundary data (Ub,0b) should be located 
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\& x (U(x),e{x))\<C „ , +C6e- cx for k = 0,1, (2.16) 



in the right region from the local stable manifold, that is, (11^,0^) should satisfy a 
condition 

U h >h s (0 h ). (2.14) 
From (12.131) . we also see that the solution z satisfies 

< c-V < z{x) < C-^— , \d k J(x)\<C—^-—r. (2.15) 
l + 5x~ KJ ~ 1 + Sx 1 x Wl (l + 5x) fc+1 v 7 

Combining ( 12. 9ft . (12. lOf) and (12.151) with using h c (z) = 0(z 2 ), we have the decay 
property of (U, 0): 

i Tsxy 

Finally we prove the existence of the solution to (12. 5p and (12.61) for the subsonic 
case M + < 1. For this case, the eigenvalues of the matrix J are Xi > and A 2 < 0, 
so that there exist a local unstable manifold and a local stable manifold. Therefore, 
the problem (12.51) and (12.61) has a solution (U, 0) satisfying (12. 7p if the boundary 
data is located on the stable manifold, that is, 

U h = h\0 h ). (2.17) 

We summarize the above observation in Lemma [2TT1 as the existence result to the 
problem ([23]) and (T23|) . 

Lemma 2.1. Suppose that \ {Uh,0\,)\ is sufficiently small. 

(i) For the supersonic case M + > 1, there exists a unique smooth solution (U,0) 
to the problem (12.51) and (12. 6p satisfying (12. 7p . 

(ii) For the transonic case M + = 1, if the boundary data (£4>,<9b) satisfies (12.141) . 
there exists a unique smooth solution (U,0) satisfying (I2.16p . 

(iii) For the subsonic case M + < 1, if the boundary data {11^,0^) satisfies (12.171) . 
there exists a unique smooth solution (U,0) satisfying (12.71) . 

The proof of Proposition 11.11 follows immediately form Lemma 12.11 Indeed, by 
using the conditions (I2.14p and (I2.17p . we precisely define the regions M° and Mr 
in Proposition [LI] as follows. Define U(u,8) and 0(u,8) by 

S:S)=-^(:-i)- < 2 ' i8) 

Note that U(x) = U(u(x),9(x)) and 0(x) = 0{u{x),B{x)) hold from (fjOJ. Then, 
defining the regions Ai° and M.~ by 

M° := {M) e M + ; U(u,6) > h s (0(u,6))}, 
Mr ■= {( U ,B) e M + ; U{u,6) = h s (0{u,6))}, 

we see that the conditions (12.141) and (I2.17P are equivalent to (11.151) and (11.171) . 
respectively. 
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2.2 Estimates for degenerate stationary solution 

The aim of the present section is to obtain more delicate estimates of the degener- 
ate stationary solution, which will be utilized in deriving a priori estimates of the 
perturbation from the degenerate stationary solution for the case M + = 1. 

Lemma 2.2. Suppose that the degenerate stationary solution exists. Namely, the 
same conditions as in Proposition W.W - (ii) are supposed to hold. Then the degenerate 
stationary solution (p, u, 9) satisfies 

05, u, 9) = (1, -1, 1) + (-1, -1, 1 - 7 )5 + 0(z 2 + 5e~ cx ), (2.20) 

{u x , ~9 X ) = 7 - l)z 2 + 0(~z 3 + 5e~ cx ), (2.21) 

\d k x {u,9)\ < C~z k+l + C5e- cx for A; = 1,2,.... (2.22) 

Proof. The estimates for (w, 9) in (I2.20p are obtained by using (12.91) . (12.101) and 




(2.23) 



which follows from (12.41) . Due to the fact that pu = —1, we have the estimate for p 
in rf2720]) . By using ([11]), fl2TT0|) and fl2TT3|) . we see that 

U x = -^-^z 2 + 0(z 3 + Se~ cx ), x = O(z 3 + Se~ cx ). (2.24) 

Differentiating ( 12.231) in x and substituting (12.241) yield the desired estimate (12.211) . 
We also have the estimates \d k (U, 0) \ = 0(z k+1 + 5e~ cx ) inductively, which give the 
estimate (12.221) due to ( 12.231) . Therefore we complete the proof. □ 



2.3 Local structure of invariant manifolds 

In order to verify the conditions (12. 14|) and (I2.17p . which ensure the existence of the 
stationary solution, it is important to make clear the local shapes of the invariant 
manifolds h c and h s . In the present section, we focus ourselves on the transonic 
case M + = 1 and show that the geometric properties of the invariant manifolds 
are characterized by the Prandtl number. In detailed arguments, we follow an idea 
in [3]. Precisely, we approximate h c and h s by polynomial functions around the 
equilibrium point as 

h%U)=c 2 U* + c 3 U 3 + 0(U*), 

h s (G) = s 2 e 2 + s 3 G 3 + 0(G 4 ) [ ' 

and obtain precise expressions of the constants Cj and Sj (i = 2, 3). 

Firstly we treat the center manifold h c . Differentiating the relation = h c (U) 
in x, we have 

O x = (h c )'(U)U x . (2.26) 

Substituting the equation (12.51) in (I2.26P and using the relation = h c {U) again, 
we have 

X 2 h c (U) + g(U, h c (U)) = (h c Y(U)f(U, h c (U)), (2.27) 
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where we have used Ai = 0. Substituting A 2 = —Cyd/(fjLK) and (12.121) in (12.271) and 
using the equalities 

e = h°{u) = o{u 2 ), (h c )'(u) = o(\u\), f(u,h c (u)) = o(u 2 ), 

we get the second order approximation of h c : 

h c (U) = ~g{U, h\U)) + 0(\U\ 3 ) = 7(T 2 ~ 2 1)2 V r " 2)£/ 2 + 0(\Uf). 

This approximation means c 2 is given by 

_ 7 ( 7 -1) 2 k 

For the case of P r = 2, that is, c 2 = 0, we compute C3 similarly as above and get 

7(7- l) 2 ^ 
C3 = — ^ > °- 

Next we obtain s 2 and S3. Differentiating U = h s (0) in x and substituting (12. 5p 
in the resultant equality, we have 

f(h s (0), 0) = (h s ) f (0) (\ 2 + g{h\0), 0)). (2.28) 

Substituting {h s )'(0) = 2s 2 + 0{0 2 ), g{h s {0),0) = O(0 2 ) and 

f(h s (0), 0) = (7 ~ ^ )V (fi - 7*)£ 2 + O(|0| 3 ), 7* := \{l 2 ' 7 + 2) > 1 
in (T2T28D . we have 

( 7 -i)V« 3 
52 = 2^ (Pr - 7 * ) ' 

If P r = 7*, that is, s 2 = 0, we also compute s 3 in the same way: 

_ 7(7 ~ l) 5 /^ 4 . n 

53 " 6d 2 > °- 

Summarizing the above observation, we have 

Lemma 2.3. Suppose that M + = 1 holds. 

(i) The local center manifold = h c (U) = c 2 U 2 + C3U 3 + 0(£/ 4 ) satisfies c 2 ^ 
and ore/?/ i/P r ^ 2. Especially, if P T = 2, i.e., c 2 = 0, the coefficient C3 is positive. 

(ii) TTie /oca/ ataftte manifold U = h s (0) = s 2 2 + s 3 <9 3 + 0(<9 4 ) satisfies s 2 | 



< 
> 

coefficient S3 is positive. 



if and only if P T > 7* := ( 7 2 — 7 + 2)/2. Especially, if P T = 7*, i.e., s 2 = 0, the 



From the local structure of the invariant manifolds in the diagonalized coordinate 
(U, 0), we obtain detailed information on the local structure of invariant manifolds 
in the original coordinate (u,8). Let 9 = h c (u) and 9 = h s (u) be a local center 
manifold and a local stable manifold in the coordinate (u,9), respectively (also see 
Figured]). Then we see that the relations 9 = h c (u) and 9 = h s (u) are equivalent to 

0{u, 9) = h c (U(u, 9)) and U(u, 9) = h s (©(«, 9)) , (2.29) 
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respectively. Therefore, substituting (12 . 1 8[) and f!2.25j) in (12.291) and solving the 
resultant equation with respect to 9, we get 

h c (u) = 1 + ( 7 - 1)(« + 1) + ~ l) (P T - 2)(« + l) 2 + 0(\u + 1| 3 ), 

2(P r + 7-I) 

h\u) = 1 - P T {u + 1) + -— ^ - 7*X W + X ) 2 + 0(\u+ 1| 3 ). 

2(P r + 7-I) 

Especially, if P r = 2 the local center manifold 9 = h c (u) satisfies 

~h\u) = 1 + (7 - l)(t* + 1) - J ( 2 ~ X \ (ti + l) 3 + 0(|n + 1| 4 ), 

i7 + 7 - 1 

while the local stable manifold 9 = h s (u) satisfies 

7(7- m ,.. , ^3 , , 1,* 



h*(u) = 1 - P r ( W + 1) + ' \ A u + 1Y + 0(\u + 1| 4 ) 

6(P r + 7 - 1) 



if P = 7*- 



3 Energy estimate 

In this section, we prove Theorem 11.21 The crucial point of the proof is a derivation 
of a priori estimates for a perturbation from the stationary solution 

(<p, ^, x ) '■= G°> - (p> s > 0)0*0 

in the Sobolev space H l . Using (jl.8p and (11.111) . we have the system of equations 

for ((p, ip, x) as 

(ft + ucp x + pV x = -{u x (f + p^), (3.1a) 

p(ll> t + M^x) + jj^2"(P _ P)x = ^xx ~ (fiU ~ pu)u x , (3.1b) 

c c ~ 1 

V PX* + -rkiP 11 ®* ~ puQx) = i^Xxx + M M x - €) - T7f(P Ux ~ PUx)- (3-lc) 



M 2 /^ x r x ' AXX ^ v x x ' M 2 + 

The initial and the boundary conditions for (ip,ip,x) follow from (II. 2p and (11.31) as 

(<P,*l>,x)(Q,x) = (<fo,ipo,Xo)(x) ■■= (Po,u ,9 )(x) - (p,u,6)(x), (3.2) 
(V,X)(*,0) = (0,0). (3.3) 

Hereafter for simplicity, we often use the notations <P := ((p,ip, x) J an d $o : = 
(<Po, ipo,Xo) J - 

To show the existence of a solution to the problem (13. ip . (13. 2p and (13.31) locally 
in time, we define a function space X(0,T), for T > 0, by 

x(o,T) := {M,x) ; v e 4 +CT/2,1+ff > (Vsx) e 4 +<7/2,2+<T > 

(p, ^ x) e C([0, T]; fT 1 (R+)), ^ G L 2 (0, T; L 2 (M + )), 

where a G (0, 1) is a constant. We summarize the existence theorem in the following 
lemma, which is proved by a standard iteration method. 
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Lemma 3.1. Suppose that the initial data satisfies 

/or a certain o G (0, 1) and compatibility conditions of order and 1. T/ien there 
exists a positive constant T , depending only on \\(Po\\b 1 +^ Q,nd \\(il>o, Xo)||b 2 + ct > such 
that the problem (13.11) . (13. 2\ and (13.31) has a unique solution (ip,il),x) G X(0,T ). 

Next we show a priori estimates of the perturbation (<£>,-?/>, x) in the space H 1 . 
Here we utilize the Poincare type inequality in the next lemma. Since this lemma 
is proved in the similar way to the paper [9], we omit the proof. 

Lemma 3.2. For functions f G /f 1 (lR + ) and w G L\(R + ) f we have 

[ \w(x)f(xY\dx<C\\w\\ Ll (f(0Y+\\f x \\ 2 ). (3.4) 

To summarize the a priori estimate, we define non-negative functions N(t) and 
D(t) by 

N{t) := sup \\$(t)\\bx, 

0<T<t 

D(t) 2 := |(^^)(t,0)| 2 + ||^(t)|| 2 + ||(^,Xx)WII^- 

Proposition 3.3. Assume that the stationary solution exists. Namely, one of the 
following three conditions is supposed to hold: (i) M + > 1 and ( 11.131) . (ii) M + = 1 
and (051) .. or (iii) M + < 1 and ffTTTj) . Let <P = {<f,ip,x) e X(0,T) &e a solution 
to (13.11) . (13.21) and (13. 3ft /or a certain constant T > 0. T/ien t/iere ezis£ positive 
constants e 3 and C independent ofT such that if N(T) + 5 < e 3 , then the solution 
$ satisfies the estimate 

\mt)\\m+ f W) 2 dr < C\\<P \\ 2 H1 . (3.5) 

We prove Proposition ^. 31 in Section I3~T1 for the case where the stationary solution 
is non-degenerate, that is, M + ^ 1. Since the decay property of the degenerate 
stationary solution for the case M + = 1 is different from that of the non-degenerate 
one, we have to modify the derivation of the estimate (I3.5P for M + = 1. It will be 
studied in Section [3721 

In deriving a priori estimates, we have to employ a mollifier with respect to time 
variable t to resolve an insufficiency of regularity of the solution obtained in Lemma 
13.11 As this argument is standard, we omit detailed computations and proceed a 
derivation of the estimates as if the solution verifies the sufficient regularity. 

3.1 Estimates for supersonic and subsonic flows 



In this section, we obtain the uniform a priori estimates of the perturbation from 
the non-degenerate stationary solution. Namely, we show (I3.5P for the case M + ^ 1. 
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In order to obtain the estimate (13.51) . we firstly derive a basic L 2 estimate. To this 
end, it is convenient to employ an energy form S defined by 

Owing to a smallness assumption on N(T), a quantity ||$||l°° is also sufficiently 
small. Hence we see that the energy form is equivalent to |$| 2 : 

c^<uj(^)<C^\ cx 2 <u;Q<C X 2 , c\<P\ 2 <£<C\<P\ 2 . (3.6) 

The solution, moreover, satisfies the follorin uniform estimates 

< c < p(t, x), 9{t, x) <C, -C < u(t, x) < -c < (3.7) 
for (t,x) G [0,T] x R + . 

Lemma 3.4. Suppose that M + ^ 1 and the same conditions as in Proposition \3.3\ 
hold. Then we have 

\mt)\\ 2 + f { V (T,0) 2 + \m x , Xx )(r)\\ 2 )dr 
Jo 

< C\\<P \\ 2 + C5 [ y x (r)\\ 2 dr. (3.8) 
Jo 

Proof. Multiplying (13.1bj) by t/>, and (13 . lc[) by xjQ, then adding up the resultant two 
equalities, we have 

{pS) t - {G? + B 1 ) x + ii-^l + kj 2 x 2 x = u x G (2) + ~e x Gf + R x , (3.9) 



G? := -(pu - pu)1> - -j^ + ~ ^X 2 , 

g? ■= j^pwty + - MTe x[(md ~ mh 

K ~ 2p _ , 

R\ ■= -jpVxXXx + -ju x xwx- 

Due to the boundary conditions fll.lOp and (13.31) . the integral of the second term on 
the left-hand side of (13.91) is estimated from below as 

-/ (G ( 1 L) + B 1 ) x dx = -(pu£)\ x=0 >ap(t,0) 2 . (3.10) 
Jr + 

In order to estimate the right-hand side of ( 13.91) . we use ( 11.141) . ( 13.41) and the fact 
|(Gi 2) ,Gf } )| < C\<P\ 2 , which follows from (EE} and Hence we have 

/ \u x G? ] + 9 X G® +R 1 \dx<C8\\(^ x , Xx )\\ 2 + C5 [ e~ cx \$\ 2 dx 
Jr + Jr + 

<CT(^(t,0) 2 +||^f). (3.11) 
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Therefore, integrating (13. 9p over (0,T) x K. + , substituting (I3.10p and (13. lip in the 
resultant equality, and then letting 5 suitably small, we obtain the desired inequality 

(GEEK. □ 

Our next aim is to get the estimate for the first order derivative (<p x , ip x , \x)- To 
do this, we first derive the estimate for (p x . 

Lemma 3.5. Suppose that M + ^ 1 and the same conditions as in Proposition ^. 31 
hold. Then we have 

\\Mt)f+ f {Mr,0) 2 + \\ Vx (r)\\ 2 )dr 
Jo 

<C\\<P \\m+C(N(t) + 5) [ D(r) 2 dr. (3.12) 

Jo 

Proof. Differentiate (13. lap in x to get 

fxt + U^Pxx + P*Pxx = /2, (3.13) 

h ■= -{2(p x if> x + 2u x (p x + 2p x ip x + u xx (f + p xx ip). 
Multiplying (I3.13P by <p x yields 

+ (~2 U ^)x = ~ P<fx ^ xx + := \ Ux{f ** + ^ 3 - 14 ^ 

On the other hand, multiplying (13.1b)) by pip x yields 

(p 2 <p x i/j)t ~ (p 2 <Ptif>)x + JpPPl = PPVx^Pxx + G 2 + R 2 \ (3.15) 

G 2 := p 3 ^. - j^p 2 f x Xx, 

R { 2 2) := -2pp x ip t 4> + p 2 ^ x (u x tp + p x i)) - —^-p(p x (9 x (p + p x x) - pu x tp x (pu - pu). 

Ml r y 

Successively multiplying (13 . 14)) by p and adding the resultant equality to (13 . 15)) . we 
have 

fv2 + P\xi) t + (fal - P 2 ^ti) x + -^fPvl = G 2 + Jfe, (3.16) 

R 2 \= pR 2 ^ + R 2 ^ . 

Owing to the outflow boundary condition on u in (11.31) . the integral of the second 
term on the left-hand side of (13. 16)) is estimated from below as 



U(fl - p 2 (ftW\ dx = -^u h ip x (t, 0) 2 > c(p x (t, 0) 2 . (3.17) 

/ X z 



Hereafter, we denote e an arbitrary positive constant and C £ a positive constant 
depending on e. The first term on the right-hand side of (13.161) is estimated as 

\G 2 \<evl + C £ \(^ X x)\ 2 . (3.18) 
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Since the second term on the right-hand side of (13.161) is estimated as 

\R 2 \ <C\%u x \<p 2 x + C5\(ip Xl rb x )\ 2 + C5e~ cx \$\ 2 1 (3.19) 
we get the estimate for the integral of R2 as 

/ \R2\dx < C(N(t) +5)D(t) 2 . (3.20) 
In deriving (13.201) . we have used the estimate 

Mfldx < \\iI)x\\l°°\\<Px\\ 2 < CWtpxWmWPxW 2 < CN(t)(\\ij; x \\ 2 H i + \\(f x \\ 2 ) 



to handle the first term on the right-hand side of ( 13.191) and the Poincare type 
inequality ( 13 Ah to estimate the third term. 

Therefore, integrating (13~T6|) over (0,T) x R + , substituting (l3T7j) . (13181) and 
( I3.20p in the resultant equality and then letting e small, we obtain 

\WA 2 + / (<Px(r, 0) 2 + ll^f) dr < C\\<P \\ 2 + C||^||^ + C f \\(i> x , X x)\\ 2 dr 
Jo Jo 

+ C(N(t)+5) [ D(r) 2 dr, 
Jo 

which yields the desired estimate (I3.12j) by substituting ( I3.8P in the second and the 
third terms on the right-hand side. These computations complete the proof. □ 

Next we estimate ip x . 

Lemma 3.6. Suppose that M + 7^ 1 and the same conditions as in Proposition ^. 31 
hold. Then we have 

ll^(t)|| 2 + f UUr)\\ 2 dT<C\\<d> \\ 2 m +C{N(t) + 5) f D(r) 2 dr. (3.21) 
Jo Jo 

Proof. Multiplying ( I3.1bl) by —ij) xx gives 

(\p^) t ~ W>M* + = + Rs, (3.22) 



1 

Mil 



G 3 := pu%b x %b xx + {9<p x + pXx)^z 



1 



-^3 := -7J?-(9xy + PxX)^xx + u x (pu - pu)%b xx - p x rb x rb t + ^-pti> 2 x . 
Notice that G3 satisfies 

\G*\ <ero 2 xx + C £ \<P x \ 2 . (3.23) 
The term R3 is estimated, by ( 11.14I) . as 

\R 3 \ <C'|^||(^,^ a; )|+C'5|(^,^)| 2 + C ( 5e- c:E |^| 2 . 
By using (13. 4p and an inequality 

/ \ii> x $M*x>il>xx)\dx< ll^lUacll^niK^,^)!! <CN(t) \\ij x \\ m ||(^,^)|| 

<CN(t)\\(<P x ,rp xx 



■ r.r l\\ ) 
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we have the estimate for the integral of R3 as 

I \R 3 \dx <C(N(t) +5)D(t) 2 . (3.24) 
Jr + 

Therefore, integrating (pj22D over (0,t) x R + and substituting (13T231) and (|3~24l) in 

the resultant equality, we obtain the desired estimate (13.211) . □ 

We finally derive the estimate for \ x . 

Lemma 3.7. Suppose that M + ^ 1 and the same conditions as in Proposition ^.^ 
hold. Then we have 

IIXx(*)f + f \\xUr)\\ 2 dr<C\\^\\ 2 m +C{N{t) + 5) [ D{rf dr. (3.25) 
Jo Jo 

Proof. Multiply ( 13. left by — Xxx to get 

(d^F PX *)t ~ (w PXxXt )x + ^ = Gi + Ra ' ^ 3 - 26 ) 



c v 1 



^j6 x (P^ + w)Xxx + ^. 
2 



R 4 := - fi(u 2 x - u 2 x )xxx + TH^x{p^ + w)Xxx + T7T7 U ^PX + 9<P)Xa 



M 2_PxXxXt+ 2M 2ptXx- 

We see that G 4 is estimated as 

\G 4 \<ex 2 xx + C £ \<P x \ 2 . (3.27) 

By a straightforward computation together with utilizing (11.141) . we see that -R4 
satisfies 

Integrating the above estimate with the aid of using inequalities 

xi Xxx 

)\dx< \\(A,Xx)\\l~\\$x\\\\{$x j Xxx) II 



(3.28) 

<CN(t)\\($ x ^ xx ,Xxx)\\\ 

I i>l\$x\ 2 dx< ||^||!oo||^,|| 2 <CiV(t) 2 ||^ll^, (3.29) 
Jr + 

we get the estimate for the integral of R4 as 

\R 4 \ dx < C(N(t) + 5)D(t) 2 . (3.30) 



Thus, integrating (13.261) over (0, t) x R + and substituting (13.271) and (13.301) in the 



resultant equality, we obtain the desired estimate (13.251) . □ 
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Proof of Proposition [3731 for M + ^ 1. Summing up the estimates (13.121) . (I3.2ip and 
(I3.25p . we have the estimate for the first order derivative <H> X as 

\\<P x (t)\\ 2 + [ {<p x (T,0) 2 + \\(<p x ^ xx , XxX )(T)\\ 2 )dT 
JO 

<C\\$ \\ 2 Hl + C(N(t)+6) [ D(t) 2 oIt. (3.31) 

Jo 

Then, adding (13. 8p to (I3.3ip and letting N(T) + 5 suitably small, we obtain the 
desired a priori estimate (13. 5p . □ 



3.2 Estimates for transonic flow 

In this section, we prove Proposition 13.31 for the case M + = 1, where the stationary 
solution is degenerate. To do this, we define a dissipative norm D(t) by 

D{t) 2 :=D(t) 2 + 6 2 [<P(t)}-2, 
where the norm [ • ] a is defined by 

[ u ] a :=([ (i + 5x) a \u(x)\ 2 dx) 1/2 . 
Using the above notation, we show the uniform a priori estimate 

f D(r) 2 dr<C\\<P \\ 2 Hl , (3.32) 



provided that N(T) +6 is sufficiently small. Since the desired estimate (13.51) immedi- 
ately follows from (13.321) . it suffices to show the estimate (13.321) . which is obtained by 
combining the estimates (13.331) and ( 13.411) . For the case M + = 1, a decay property 
of the degenerate stationary solution is worse than the non-degenerate stationary 
solution. Therefore, in deriving L 2 estimate of $ summarized in Lemma 13.81 we 
have to utilize the precise estimate (12.211) of the degenerate stationary solution in 
order to estimate the term u x G\ + 9 X G\ in (13.91) . 

Lemma 3.8. Suppose that M + = 1 and the same conditions as in Proposition 13.31 
hold. Then we have 

\mt)\\ 2 + f (<p(r, 0) 2 + 5 2 [<P(t)} 2 _ 2 + ||(^,x,)(r)ir) dr 
Jo 

<C\\<P \\ 2 + C5 [ y x (r)\\ 2 dT. (3.33) 
Jo 

Proof. Notice that the solution (p, u, 9) satisfies 

(p,u,9) = (1,-1,1) + 0(N(t)+ 6), (3.34) 

which follows from fl2T20l) and ||^(t)||L<» < CN(t). Using the property (13T34D and 
(I2.20p . we see that G\ is divided into a main quadratic form and residue terms as 

G? = -i? - -x 2 + - ~<px + 0(N(t) + 6)\<P\ 2 . 

Ill 
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Hence we see from the above expression and (12.211) that the first term on the right- 
hand side of (I3.9P satisfies 

t(2) _ 7 + 1-2/ ,,2 , 1, 2 7-1 , , 1 



w^i tttZ [ifj +- X <P$ + ~<PX) 

2d V 7 7 7 / 

+ 0(iV(t) + 5)~z 2 \$\ 2 + 0(5)e- ca: |^| 2 . (3.35) 



By a similar computation, we have 



Gf = --^ 2 + |x 2 + - c v ^X + 0(N(t) + 5)\<P\ 2 , 
where we have also used the fact that 

u(s)= 1 -(s-lf + 0(\s-lf). (3.36) 

Therefore, due to ( 12.211) . the second term on the right-hand side of ( 13.9f) satisfies 
S 7 + l~ 2 /7-l 2 1 2 1 , 1 , 

+ 0{N{t) + 5)~z 2 \$\ 2 + 0{5)e- cx \<P\ 2 . (3.37) 
Summing up the expressions (13.351) and (13.371) . we have 

u x G® + ~6 x Gf = -^i 2 ^, ^ x ) + 0(N(t) + 5)~z 2 \<P\ 2 + 0(5)e~ cx \<P\ 2 , 

47a 

(3.38) 

F x (<p, V, X) ■= (7 - IV 2 + 2 7 ^ 2 + X 2 ~ 2(7 - 1)W> + 2^X- 
The quadratic form Fi(<p, ip, x) is positive definite since 

F 1 (<p, ^ X ) = (7 - 1)(V - ^) 2 + (V> + X? + 7^ 2 > c|^| 2 . (3.39) 
Due to (I2.2ip . the remaining term Ri is estimated as 

\Ri\ <C6(z 2 \<P\ 2 + \(^ Xx )\ 2 + e- cx m 2 ). (3.40) 

Therefore, integrating flU} over (0,t) x R + , substituting (ET5gjl . (13391) and dSHOD 
in the resultant equality and letting iV(£) + 5 suitably small, we obtain 

||<2>|| 2 + [\<p( T , 0) 2 + 5 2 [<P] 2 _ 2 + ||(^, XxOf) cfr 
Jo 

<C||^oll 2 + ^/ [ e- cx \$\ 2 dxdr, 



where we have used ( 12.151) . Finally, to estimate the last term on the right-hand side 
of the above inequality, we utilize the Poincare type inequality (13.41) . Consequently, 
we arrive at the desired estimate f!3.33|) and complete the proof. □ 

Next we show the estimate for the first order derivative <& T . 
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Lemma 3.9. Suppose that M + = 1 and the same conditions as in Proposition 13.31 
hold. Then we have 

ll<M*)f + / {f x (r,0) 2 + \\{<f x ,^ xx ,Xx X ){r)\\ 2 )dT 



<C\\<P \\ 2 H i+C{N{t) +5) [ D{rfdr. (3.41) 

Jo 



Proof. In the present proof, we only show the estimate for the remained terms R2, 
i?3 and R4. The other part of the derivation of (I3.4ip is almost same as that of the 
non-degenerate case, so we omit the details. Using ( 12. 22ft . we see 

1(^,^3,^)1 <c\(i> x , Xx mM^,^,x xx )\ + c^ x \ 2 

+ C5~z 2 \$\ 2 + CS\(<P X , ?p xx , xxx)\ 2 + C5e~ cx \<P\ 2 . 
By computations similar to ( 13.281) and (I3.29f) . and by using ( 12.151) . we have 

1(^2,^3,^4)1 dx < C(N(t) + 5)D(t) 2 . 



Therefore, following the same procedure of the derivation of (13.121) . (13.211) and (13. 25j) 
and using the above estimate for the remaining terms, we obtain the desired estimate 
(EH}. □ 



3.3 Proof of Theorem O 

This section is devoted to the proof of Theorem 11.21 Firstly we prove the existence 
of the solution in the sense of (I1.18P globally in time. Since the existence time T Q in 
Lemma 13.11 depends on the Holder norm of the initial data, we have to show the a 
priori estimate in the Holder norm. Precisely we prove 

||p|| s i+ CT /2,i+<T, \\(u, 9)\\ B i+*/2,2+v < C{T) (3.42) 

for the solution (p,u,9) satisfying (p,u,9) — (p,u,9) G X(0,T), where C{T) is a 
positive constant depending on T, ||po||s 1 + CT , ||( M o, #o)||# 2 + ct an d IKv^o, V'o, Xo)!!^ 1 - To 
obtain the estimate (I3.42p . we rewrite the system (jl.8p in the Eulerian coordinate 
into that in the Lagrangian mass coordinate, and then apply the Schauder theory 
for parabolic equations studied in [4] with the aid of the H l uniform estimate (13.51) . 
Since the derivation of the Holder estimate (I3.42p is same as that in [H] studying 
the stability of the stationary solution for an isentropic model, we omit the details 
of the proof. Therefore, combining Lemma 13.11 and the estimate ( 13.421) by using the 
standard continuation argument, we obtain the existence of the solution globally in 
time. Moreover, we see that the solution verifies 

sup pK*)|&i+ / D{t)dt<C\\$ Q f m . (3.43) 

te[o,oo) Jo 

Next we show the stability (I1.19p . For this purpose, it suffices to show that 

\\$ x {t)\\ ^0 as t -> 00 
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since we see ||^(t)|| £00 < C 1 1 ^ C*) 1 1 1/2 1 1 (*) 1 1 1/2 and \\$(t)\\ < C due to the H 1 
uniform estimate (I3.5p . Let I(t) := \\(p x (t)\\ 2 . By a similar computation to [9], we 
have 



I(t) < CD(tf 



which gives 4/ G -^ 1 (0, oo) owing to (13.431) . Combining this fact with / G -^ 1 (0, oo), 
which is a direct consequence of (I3.43p . we have I(t) — > 0, i.e., H^COII as 
t — > oo. The convergence IKV'auXxX^H ^ is proved in the similar computations. 
Consequently, we prove (11.191) and complete the proof of Theorem 11.21 



4 Weighted energy estimate 

In this section, we show the proof of Theorem 11.31 Precisely, we obtain conver- 
gence rates of the solution toward the stationary solution by using a time and space 
weighted energy method. 



4.1 Estimates for supersonic flow 

This section is devoted to showing the convergence ( 11.201) for the case M + > 1. To 
this end, we define weighted norm E a (t) and D a (t) by 

E a (t) 2 := \Mt)\\h + W)\l 
D a (t) 2 :=D(t) 2 + a\<P(t)\l-i + \(^X x )(t)\ 2 a 

and obtain the weighted energy estimates summarized in the next proposition. 

Proposition 4.1. We assume that M + > 1 and ( 1 1 . 1 3 j) hold. Let <P = (</?, Vsx) G 
X(0,T) be a solution to (JET}, (Q and (J33J) satisfying <P G C([0, T]\ L 2 a {R + )) for 
certain constants a > and T > 0. Then there exist positive constant £4 and 
C independent of T such that if N(T) + 5 < £4, then the solution <& satisfies the 
following estimates 

(1+tyE^it) 2 + f (l + rYD^irfdr < CE a (0) 2 , (4.1) 
Jo 

for an arbitrary integer j = 0, . . . , [a] and 

(1 + t)*Ea{t) 2 + f\l + rfD {r) 2 dr < CE a (0) 2 (l + t)^ a (4.2) 
Jo 

for an arbitrary £ > a. 

The proof of Proposition 14.11 is based on the time and space weighted estimate 
of <P in L 2 (R + ) and the time weighted estimate of <P X . 
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Lemma 4.2. Suppose that the same conditions as in Proposition \4:.l\ hold. Then we 
have 

(i + + Ai + far, o) 2 + p\*{r)\i-i + \y>., x-)wi?) & 

Jo 

<C\<P \l + CZ r(l + r)^Xr)||dr + C5 / (1 + t)% x (t)\^ dr (4.3) 
Jo Jo 

for arbitrary (3 G [0, a] and £ > 0. 

Proof. Multiplying ( 13. 9ft by a weight function w(t,x) := (1 + t)^(l + x)' 3 , we have 



(u;p£) t - {w(G? + B 1 )} x + w x G? + w^-^ 2 x + 



# 2 ' 

= w tP S - w x B 1 + w{u x G ( ? ] + 9 x Gf ] + Ri). (4.4) 

The integral of the second term on the left-hand side of (14.41) is estimated from below 

as 

-/ {wiG^+Bt^Jx^cil+tf^O) 2 . (4.5) 

Due to (12.201) . ( 13.341) and (13.36}) . the term G^ is divided into a quadratic form and 
remaining terms as 

G? = 2M 2 1 \ J _ 1) F ^ 4>> X) + 0(N(t) + 5) \<P\ 2 , (4.6) 
F 2 (<p, V>, x) ■= (7 - l)^ 2 + M 2 7 ( 7 - l)^ 2 + x 2 - 2( 7 - 1)(^ + x)V>. (4.7) 

Notice that the quadratic form F2 is positive definite owing to the assumption M + > 
1 since 

F a (^V,x) = (7-l)(^-^) 2 + {(7-l)^-x} 2 +7(7-l)(M 2 -l)^ 2 > c\$\ 2 . (4.8) 

Thus, substituting (14.81) in (I4.6p . we have the estimate of the third term on the 
left-hand side of (14.41) from below as 



I w x GP dx>{c- C(N{t) + 5)}P(1 + t) € \$\}_ v 
Jr+ 



(4.9) 



The first and second terms on the right-hand side of (14.41) are estimated with the 
aid of the Schwarz inequality as 

/ \wp£\dx<Ci{l+tf- l \<P\% (4.10) 

[ \wMdx <C(3{\ + ^{6^1^ + C £ \{^ Xx )\l-i)- (4.H) 
Jr + 

In the similar way to the derivation of ( 13. lip , we estimate the remaining terms in 
ItPj) . by using (fTTl) and (E3D, as 

I w\u x G {2) + 9 x Gf ] + R 1 \dx<C5(l+tY(v(t,0) 2 +\\<P x \\ 2 ). (4.12) 
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Therefore, integrating (fl~4"|) over (0,t) x R + , substituting f)4~5l) and (J4~9l) - (14.121) in 
the resultant equality and letting e and N(t) + 5 sufficiently small, we arrive at 

(1 + tf\$% + l\l + T )*(<p{T, 0) 2 + mu + kv>«, dT 

Jo 

<c\<P \l + ct [ (l+rf-^dT+cs [ (i + T)ty x \\ 2 dr 

Jo Jo 

+ C(3 [ {l+rf\{^x,)\l-xdr. 
Jo 

We finally apply induction with respect to (3 to estimate the last term on the right- 
hand side of the above inequality. This computation yields the desired estimate 
(14. 3D . Consequently, we complete the proof. □ 

Letting (3 = in (14. 3p . we have the time weighted estimate 

(1 + t)«||*(t)|| a + Al + r)^(r, Of + \\(^ Xx )(r)\\ 2 ) dr 
Jo 

<C\\4> \\ 2 + Ct f(l + Tf- l \\3>{T)\\ 2 dT + C5 [\l + T^y x (r)\\ 2 dT (4.13) 
Jo Jo 

for an arbitrary £ > 0. 

We state below the time weighted estimate for the first order derivative <P X . Since 
the proof of this estimate is almost same as that of (I3.3ip . we omit the details and 
only summarize the result in the next lemma. 

Lemma 4.3. Suppose that the same conditions as in Proposition \A.l\ hold. Then we 
have 

(1 +^||^.(t)|| 2 + Al + r)« (^(r,0) 2 + ||(^,^,X,,)(r)|| 2 ) dr 
Jo 

<C||^ || 2 J1 +^ [\l + rf-'imr)]^ dr 
Jo 

+ C(N(t) + S) [ {l + TfD{rfdr (4.14) 
Jo 

for an arbitrary £ > 0. 

We conclude this section by giving the proofs of Proposition 14.11 and Theorem 

ra-(i). 

Proofs of Proposition 14.11 and Theorem 11.31 - (i) . Summing up the inequalities (14. 13|) 
and (I4.14p . we have the time weighted H l estimate 

Jo 

<Cp> ||^ + C£ !\l+tf-^{r)f H ,dr. 
Jo 
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Add (14. 3 p to the above inequality to obtain 

(1 + tfE,{tf + f\l + rf{mr)\U + D,{rf) dr 
Jo 

< CEpQSf + C£ [\l + r^d^r)!? + D,{rf) dr, 
Jo 

where we have used the inequalities 

By applying induction with respect to ft and £, studied by [8] and [16], we have 
the desired estimates (14.11) and (14.21) . The convergence ( 11.201) immediately follows 
from (14. 2p and the Sobolev inequality. Consequently, we complete the proofs of 
Proposition 14.11 and Theorem 11.31 - (i). □ 

4.2 Estimates for transonic flow 

In this section, we show the convergence (ll.2ip for the case M + = 1 by deriving the 
time and space weighted estimate in H l . To do this, we define weighted norms by 

N a (t) := sup E a (r), E a (t) := |[#(t)]|i, a , 

0<T<t 

where |[-]| S)Q is the s-th order Sobolev norm corresponding to [•]„: 

i[«ii.,«:= (DM) 1/2 = (E / o-+8zr\dZu(x)\*dx) 1/2 . 

k=0 k=0 

Proposition 4.4. We assume that M+ = 1 and (fl~T5|) hold. Let $ e X(0,T) 
be a solution to (EHJ), (Q and (ED satisfying <P E C([0, T]; H^(R + )) for certain 
constants a G [1, 2(1 + y/2)) and T > 0. ITien i/iere ea;is£ positive constants £4 and 
C independent ofT such that if 5~ 1 / 2 A r Cf (T) + 5 < £4, t/ien t/ie solution $ satisfies 
the following estimates for t € [0, T] : 

(1 + t) j E a ^(t) 2 + [ (1 + ryb^rfdr < C<T 2i £ a (0) 2 (4.15) 

/or an arbitrary integer j = 0, . . . , [a/2] and 

(1 + t)^o(t) 2 + f (1 + r) f A)(r) 2 dr < C5- a E a (0) 2 (l + t)^ 2 (4.16) 

for an arbitrary £ > a/2. 



In order to prove Proposition 14.41 we have to derive time and space weighted 
estimates not only for <& in L 2 but also for the first order derivative <& x . In deriving 
the weighted estimate for <P in L 2 , we utilize the following interpolation inequality 
to handle several nonlinear terms. 
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Lemma 4.5. Let (3 > 1. T/ien a function f G i7g(R+) satisfies 

(1 + fe)^ 1 !/^)! 3 ^ < ^- 3/2 [/]i(/(0) 2 + 5 2 [f} 2 _ 2 + [/.]?)■ (4.17) 



Since we can prove ( 14.171) in the similar way to the paper [19J, we omit the proof 
of Lemma [4.51 For details, see Lemma 5.1 with p = 2 and a = (3 in [19|. 

Then we show the time and space weighted L 2 estimate. In deriving this estimate, 
we have to assume that the weight exponent a is less than 2(1 + y/2) in order to 
obtain the dissipative term S 2 [^>} 2 ^_ 2 . Moreover, to control the term 2r +1 |<2>| 3 in 
f )4.2ip . we have to assume the smallness of Hence we need a condition a > 1, 
too. 

Lemma 4.6. Suppose that the same conditions as in Proposition ^ .41 hold. Then we 
have 



Jo 

<Cl<P ] 2 8 + CZ f\l + Tf- 1 [<P{r)]ldr^C{5- 1 l 2 N p {t)+5) f\l + TfD p {r) 2 dT 
Jo Jo 

(4.18) 

for arbitrary constants (3 G [1, a] and £ > 0. 



Proof. In the present proof, we employ a spatial weight function 

Ad \V2 
,7 + 1, 

Notice that iu ~ <J~0(1 + 5x)' 3 holds due to fl2TToD . Multiplying <^Mj by the weight 
function w(x), we have 



w(x) := h 2 z(x)-^ h:=(——\ 

W + 1/ 



(pwS) t - {w(G? + B 1 )} x + w x G? + w a B x + w(fJf% + k^ x 2 x ) 

= w{u x G {2) + 6 x Gf) + wRl (4.19) 

The remainder of the present proof is divided into three steps. 
Step 1. Firstly we show that the equality (14.191) is rewritten as 

(pwS) t - {w(G? + B 1 )} x + F = R, (4.20) 

where F is defined by 

F := jW x F 2 + ~z-? +2 F 3 + ii (hip x + * 2-* + k (hx* + j*x) V", 

F 3 = F 3 (<p, X ) := -Ftfr, V, X) + -Fi{<f, V>, X) ~ j^(^ 2 + ^), 
7 7 ft/ 

F 4 = F 4 (</?, V, X) := (3 - 7^ 2 + X 2 + 2( 7 - 1) W> + 2</> X 



Stationary wave to viscous heat-conductive gases 



29 



and the remaining term R satisfies 

\R\ < c(Np(t) + 5)(rP +2 \<P\ 2 + rP\^ x , Xx )\ 2 )+C5e- cx rP\<Pf 

+ Cz' (3+1 \<P\ 3 . (4.21) 

For this purpose, we show that the third term on the left-hand side of (14.191) verifies 
a decomposition 

Wx G<p = ^w x F 2 + -pz~P +2 F A + 0(\<P\ + z 2 + 5e- cx )z-P +1 \$\ 2 , (4.22) 
2 7 

where F 2 is defined in (14.71) . Using the fact that 

(p, u, 9) = (1, -1, 1) + (-1, -1, 1 - 7 )5 + 0(\<P\ + ~z 2 + 5e~ cx ), (4.23) 
which follows from (I2.20p . we see that the terms in satisfy 

-pu£ = ^ 2 + \p 2 + |x 2 + (^V + ^X 2 )~z + 0{\*\ +z 2 + 5e-m 2 , 

-(p-p)V = --<{*!) - ~^X+ (~ — -^0 + ~^x)z + 0(\$\ +z 2 + 5e- cx )\$\ 2 . 
7 7 ^7 7 ' 

Summing up the above two equalities, we see that G^ satisfies 

= ^F 2 + ^zF A + 0(\$\ + z 2 + 5e~ cx )\<P\ 2 . (4.24) 
2 27 

Furthermore, by differentiating the weight function w(x) and using (12.131) . we have 

w x = 2(3z~ fi+l + 0(/3^+ 2 ). (4.25) 

Multiplying (I4.24p by (I4.25P yields the desired equality (14.221) . We also see that the 
fourth and the fifth terms on the left-hand side of ( 14. 19ft are rewritten as 

( Q ,2 2\ 

w x B t + w [VqVs + k-^Xx) 

+ o(\$\ + ~ z ){r? +2 \<P\ 2 + r'KVfe, x,)\ 2 ), (4.26) 

which follows from (12.201) . (I4.23P and (I4.25p . To estimate the right-hand side of 
flUHD, we use (jSSHD and (jS3QD to obtain 

w(u x G? + 9 X G?) = - -r^Frifp^x) + 0(N(t) + 5)r^ 2 \<P\ 2 

7 (4.27) 

+ 0{5)e- cx z~ l3 \<P\ 2 , 
\wRil < C5(z~P +2 \<P\ 2 + z~P\{tfj x , Xx )\ 2 + e- cx z~P\<P\ 2 ). (4.28) 

Therefore, substituting (14T22D . (14T26]) . (14371) and (jQBD in (|Q5jl . we obtain the 
desired equality (14.201) . 

Step 2. Our next aim is to show that F satisfies the estimate from below as 

F>cS-^ 2 \<P\ 2 + c^\(ij x , Xx )\ 2 (4.29) 
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provided that (3 G [0,2(1 + y/2)). Let A 2 be a real symmetric matrix satisfying 
F 2 = $TA 2 <P, i.e., 

/ 7 -l I-7 \ 
A 2 := 1-7 7(7-1) I-7 . 
V I-7 1 / 

We see that the matrix A 2 admits three distinct eigenvalues 0, v_ and u + satisfying 

v± = i(7 2 ± a/V _ 4 7 3 _|_ i2 7 2 _ 2O7 + 12) and < i/_ < 

Let gi, ^2 and 53 be unit eigenvectors of A 2 corresponding to the eigenvalues 0, i/_ 
and u + , respectively. Especially, we obtain 

gi = (1, 1,7- l) T q where q := {~f 2 - 2 7 + 3)~ 1/2 . 

Furthermore, we employ a new function $ defined by 

where Q := (91,52,93) is an orthogonal matrix. Using the fact that Q J A 2 Q = 
Q~ 1 A 2 Q = diag(0, z/i, z/ 2 ), we see that the quadratic form F 2 satisfies the estimate 
from below as 

F 2 = (QS) t A 2 QS = uj 2 + v +X 2 > c|(^,x)| 2 . 

Combining this estimate with the inequality w x > c(3z~^ +1 , which follows from 
P5| with 8 < 1, we have 

|^F 2 >c^ +1 |(^,x)| 2 . (4.30) 

Next we employ a real symmetric matrix A 3 satisfying F 3 = <P J A 3 <P. Let A 3 := 
(a-ij)ij '■= Q T A 3 Q. Then we see that 

F 3 = (QS) T A 3 QS = £tA 3 $ = a u 2 + 0(|(^, X )\ 2 + \<f>(4> + X)D- (4-31) 
Since the sign of an will play an important role later, we obtain it explicitly: 

fin = <Zi"^3<Zi = ^s|*=« = ^9 2 (4 + 4/5 - /3 2 ). (4.32) 

Owing to the above observations, we show that the first and the second terms in 
the definition of F satisfy 

jW x F 2 + z~P +2 F 3 > cz- p+2 \$\ 2 (4.33) 

provided that (3 G [0, 2(1 + a/2)). Notice that the estimate (I4.33P immediately yields 
the desired estimate (I4.29p . If (3 — 0, the quadratic form F 3 is positive definite, i.e., 
F3 > c\<D\ 2 since we have F 3 = -F1/7 and the positivity of F\ due to (I3.39p . Thus, 
owing to the continuous dependency on /3, there exists a positive constant /3* such 
that F 3 > c\<P\ 2 holds for (3 G [0, /?*], where c is independent of /3. Namely, (I4.33P 
holds for (3 G [0, /?*]. 



Stationary wave to viscous heat-conductive gases 



31 



Next, we show f!4.33j) for (3 G [{3*, 2(1 + v2)). Note that the constant an is 
positive due to (14.321) . Thus, using (I4.30P and (I4.3ip . we have 

^w x F 2 + z^ +2 F 3 

> c(3^ +1 \^,x)\ 2 + a u r^ 2 2 - Cr^ 2 {\(ij,x)\ 2 + \<p$ + x)\) 

which yields (14.331) if 5 is sufficiently small. Therefore, we have shown that the 
estimate (I4T33D holds for (3 G [0, 2(1 + v 7 ^))- 

Step 3 . Finally we prove (T4~T8j) by using (T4~20l and (T4~29i) . Using (l2~T5j) . we have 
the estimate for the integral of the second term on the left-hand side of (14.201) as 

-/ {w(G? + B 1 )} x dx>cS-^(t,0) 2 . (4.34) 

Furthermore, due to the Poincare type inequality (13.41) and the inequality (14.171) . 
integrating (I4.2ip yields 

f \R\ dx < C5-P{5- 1/2 N p (t) + 5)D p (t) 2 , (4.35) 

where we have used [$(£)] i < Np(t) for (3 > 1. Consequently, integrating (I4.20p 
and substituting (I4.29p . (14.341) and (14.351) in the resultant equality gives the desired 
estimate (14.181) . Thus we complete the proof. □ 

Next we show the estimate for the first order derivative <P X . Owing to the de- 
generate property of the transonic flow, we have to employ the spatially weighted 
energy method for the estimate for (p x . 

Lemma 4.7. Suppose that the same conditions as in Proposition ^ .41 hold. Then we 
have 

(l + t)Z[<P x (t)}l+ [ (l + r)«(^(r,0) 2 +[^.(r)] 2 + [(^, Xra )(r)] 2 )rfr 
Jo 

<C\m\ p + Ci f\l+T)^\[<l>{T)}\l,dT 
JO 

+ C(5- 1 / 2 N (3 (t)+6) [ (l + rfD^rfdr (4.36) 
Jo 

for arbitrary constants (3 G [1, a] and £ > 0. 

Proof. We only show the estimate for ip x as the other estimates for (ip x , Xx) can be 
established by similar computations. Multiplying (I3.16P by a spatial weight function 

w := (1 + dx) 13 , we get 

{ w [^<fl + (?<PxV) } + \w - P 2 <Pti>) } + wp<f 2 x 

= wG 2 + wR 2 + w x ~ P 2 VtV) ■ (4.37) 
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The last term on the right-hand side of (14.371) is estimated as 

dx < e[<p x ]% + C £ (5 2 [$}y 2 + + C5Dp{t)\ (4.38) 



w. 



P 2 2 / 

2 U( fx~ P <PtV> 



where e > is an arbitrary constant and C e is a positive constant depending on e. 
The other terms in (14.371) are estimated in a same way as the proof of Lemma 13.51 
For instance, the remaining term R% verifies the estimate 

/ w\(R 2 ,R 3 ,R 4 )\dx<C(N p (t)+6)ba(t) 2 , (4.39) 
which follows from the Poincare type inequality (13.41) and the inequality 

xxi Xxx )| dx < \\(lpx,Xx)\\L°°l$x}f3[($x,?Pxx,Xxx)}t3 
<CNp(t)[^x,^xx,Xxx)}l. 

Therefore, integrating (071) over (0,t) x R +> substituting fl3TT8|) . (OSI) and (1439]) 
in the resultant and then letting e sufficiently small with using (I4.18p . we obtain the 
estimate for ip x in (14.361) . The estimates for {ip x ,Xx) are obtained by similar com- 
putations to Lemma [3761 and [3~71 with using the estimate (I4.39p . Thus we complete 
the proof of the desired estimate (I4.36p . □ 

Proofs of Proposition 14.41 and Theorem 11.31 - (ii) . Summing up the estimates (14. 18j) 
and (I4.36p . and letting S^^Npit) + 5 be suitably small, we have 

(1 + tf\m)]\l, + f (1 + rf{5 2 [<P{T)\l_ 2 + D p {rf) dr 



o 



< C\ml, + d /'(I + rf-\[<P{r% + D,{rf) dr, 
Jo 

which yields the desired estimates (14.151) and (14.161) by an induction with respect to 
P and ^ (see [IS]). The convergence rate (ll.2ip follows from the estimate (I4.16P with 
the aid of the estimate < |[^]|f Q < ||^||^i- We consequently complete the 

proofs. □ 
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